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A fully nonlinear, time-asymptotic theory of resonant par- 
ticle trapping in large-amplitude quasi-parallel Alfven waves 
is presented. The effect of trapped particles on the nonlin- 
ear dynamics of quasi-stationary Alfvenic discontinuities and 
coherent Alfven waves is highly non-trivial and forces to a sig- 
nificant departure of the theory from the conventional DNLS 
and KNLS equation models. The virial theorem is used to 
determine the time-asymptotic distribution function. 
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The magnetic fluctuations frequently observed in So- 
lar Wind and Interstellar Medium plasma have been the 
subject of protracted and intense observational and the- 
oretical scrutiny. It is likely that these fluctuations are 
nonlinear Alfven waves, in which the ponderomotive cou- 
pling of Alfvenic magnetic field energy to ion-acoustic 
quasi- modes has modulated the phase velocity va , and so 
caused steepening and formation of discontinuities . 
Such rotational and directional discontinuities have in- 
deed been observed in the Solar Wind, and are proba- 
bly quasi-stationary waveform remnants of nonlinearly 
evolved Alfven waves 0] . 

Beginning with the work of Cohen and Kulsrud 
the theory of quasi-parallel, nonlinear Alfven waves has 
received a great deal of attention and has spawned 
in a variety of modifications of the wave envelope evo- 
lution equation, referred to as the Derivative Nonlinear 
Schrodinger (DNLS) equation. However, almost all at- 
tention has been concentrated on developing and extend- 
ing the fluid theory of such waves, leaving issues of par- 
ticle kinetics aside. Nevertheless, some attempts to in- 
corporate particle dynamics into the DNLS model have 
been made, both analytically (7| and (very extensively) 
via particle- and hybrid-code simulations (§] . Progress in 
constructing an analytical kinetic-MHD model of non- 
linear coherent Alfven waves occured recently by the 
self-consistent inclusion of linear Landau damping [^]|| 
and gyro-kinetic (e.g., ion-cyclotron) effects B. How- 
ever, even in these treatments, wave-particle resonant in- 
teraction is treated perturbatively and calculated using 
the linear particle propagator. This technique fails for 
a large-amplitude wave propagating in a finite-/? plasma 
(here f3 is the ratio of kinetic and magnetic pressure) 
because of non-perturbative effects associated with par- 
ticle trapping in the field of the wave. In this Letter, 



we extend the theory of 'kinetic' nonlinear Alfven waves 
to the strongly nonlinear regime where trapped particles 
are important. 

In the finite-/3, isothermal regimes typical of the Solar 
Wind (i.e., c s ~ va, T e ~ Tj) at 1 AU, resonant inter- 
action of the plasma with ion-acoustic quasi-modes is a 
critical constituent of the wave dynamics. The very exis- 
tence of rotational discontinuities is due to the nonlinear 
coupling of Alfven waves to (linear) Landau dissipation 
Q. Here, linear Landau dissipation refers to damping 
calculated perturbatively, assuming a Maxwellian parti- 
cle distribution function (PDF), and thus with a time- 
independent rate coefficient. This mechanism enters the 
Alfven wave dynamics nonlinearly (i.e., in proportional to 
the magnetic energy density of the wave train) because it 
enters a functional with the parallel ponderomotive force 
cx d z (Bj_/8ir). The 'kinetic' wave equation, called the 
Kinetic Nonlinear Schrodinger (KNLS) equation, is 



v"a d 2 b 



|b(z')| 2 dz' 



(1) 



where b = (B x + iB y )/Bo is the normalized complex 
wave amplitude, f2 c is the ion-cyclotron frequency, the 
coefficients mi and m,2 are functions of (5 and T e /Ti only 
(see H), and V means the principal value integration. 

Obviously, particles which are near resonance with 
the wave [v ~ va) will be trapped by the ponderomo- 
tive potential (or equivalently, by the electrostatic fields 
of driven ion-acoustic perturbations). Particle bounce 
motion significantly modifies the PDF near resonance, 
since trapped particle phase mixing results in flattening 
of the PDF (for resonant velocities) and formation of 
a plateau. Thus, the linear calculation of the Landau 
dissipation, while correct for times short compared to 
the typical bounce (trapping) time, r <C r tr , fails for 
quasi-stationary waveforms for times r > tjvx 3> r tr 
(t~nl is the typical nonlinear wave profile evolution time) . 
Hence, Landau dissipation should be calculated non- 
perturbatively to determine the resonant particle response 
to the nonlinear wave. 

Of course, the nonlinear Landau damping problem is, 
in general, not analytically tractable, as it requires ex- 
plicit expressions for all particle trajectories as a function 
of initial position and time. Such trajectories cannot be 
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explicitly calculated for a potential of arbitrary shape. 
Usually, a full particle simulation is required to obtain 
this information. In some cases, an approximate ana- 
lytic expression for the wave profile shape is known and 
may be assumed to persist, while the wave amplitude 
varies. Calculations defined in this way has been imple- 
mented for the special cases of sinusoidal P] and solitonic 
JlO| wave modulations. Other approaches either seek the 
asymptotic (r — *■ oo) PDF for a given (undamped) wave- 
form or exploit the universality of the process of 
de-trapping of resonant particles from a wave potential 
of decreasing amplitude OJl3| |. These approaches, how- 
ever, do not appear to be useful for the problem consid- 
ered here. 

The goal of this work is to investigate how trapped 
particles modify nonlinear wave evolution, assuming no 
restrictions on the shape of the wave-packet modulation. 
Thus, the motion of particles is treated self-consistcntly. 
We show that, in the two important limits of short-time 
(t -C Ttr) and long-time (r 3> r tr ) evolution, the prob- 
lem admits analytic solutions. In the limit r <C r tr , we 
recover conventional linear Landau damping. This sup- 
ports the validity of the KNLS theory as a means for 
studying the emergence of Alfvenic discontinuities. In 
the opposite limit r 3> Ttr, the virial theorem is used 
for determination of the time-asymptotic trapped parti- 
cle response. Although the damping rate vanishes due to 
phase mixing, the effects of trapped particles are highly 
non-trivial, leading to a significant departure of the the- 
ory from the familiar form of the DNLS and KNLS mod- 
els. First, the power of the KNLS nonlinearity associated 
with resonant particles increases to fourth order when 
trapped particles are accounted for. Second, the effective 
coupling now is proportional to the curvature of the PDF 
at resonant velocity, /o'(ua), and not its slope, /o(va), as 
in linear theory. Third, the phase density of trapped par- 
ticles is controlled by the plasma f3. Finally, we combine 
these to obtain the wave evolution equation which gov- 
erns the long-time dynamics of quasi-stationary Alfvenic 
discontinuities. The equation is the principal result of 
this Letter. 

We should state here that particle trapping may be ab- 
sent in higher than one dimension. Indeed, for k^pi 3> I 
is the perpendicular component of the wave vector 
and pi is the ion Larmor radius), then the longitudinal 
Cerenkov resonance to = knvu is satisfied for all parti- 
cles having v = va/ cosG, but with v±_ arbitrarily large. 
Thus, all particles with velocities v > va interact with 
a wave and a plateau cannot form, while a non-thermal 
tail of energetic particles may result instead. However, 
if the magnitude of the ambient magnetic field is strong 
enough so that k±pi <C 1, quasi-one-dimensionality is re- 
covered. This last situation is, in fact, typical for waves 
propagating in the Solar Wind. 

For reasons of notational economy, let's introduce the 
trapping potential U(z) = B^/8irno, where no is the 



unperturbed plasma particle density. Then, the char- 
acteristic bounce frequency [|| in our case is t^} ~ 
ky/U/rrii {rrii is the ion mass). The characteristic non- 
linear frequency at which the wave profile changes ap- 
preciably is readily estimated from Eq. (0) to be ~ 
m\kvA(Ji\i 'By). From comparison of these two time- 
scales, we conclude that the wave potential, as seen by 
a trapped particle, is steady-state (i.e., roughly constant 
on the particle bounce time) when t^l 3> Ttr, so that 



B ± /B 



< 



1. 



(2) 



That is, particle phase mixing is very efficient for weakly 
nonlinear waves. Note, this condition (||) is consistent 
with the derivation of the KNLS, for which B±/Bq <C 1 
(weak nonlinearity) is assumed. Let's now rewrite Eq. 
([0) in a generic form: 
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Here Sn^R is the density perturbation due to the non- 
resonant (bulk) response of the PDF. It is roughly pro- 
portional to |b| 2 . Shr is the resonant particle contribu- 
tion. It is responsible for strongly nonlinear feedback via 
the distortion of the PDF by a wave. It was also re- 
sponsible for linear damping in the KNLS equation. It is 
interesting that the very possibility to write the general- 
ized KNLS equation in the form ([}]) relies on the intrinsic 
time reversibility of the Vlasov equation, linear or nonlin- 
ear. Indeed, one can formally write the resonant particle 
response as 



5n R (x X||£ [U(z)} 



(4) 



where K. is some normalized kinetic operator acting on a 
wave field. Time reversibility implies K-K- = — I, see ||. 
This fact has been crucial for the derivation of KNLS. 
The constant xil plays a role of effective dissipation coef- 
ficient (thermal conductivity) in the linear Landau damp- 
ing theory. 

The resonant particle response is calculated using Li- 
ouville's theorem, which states that the local PDF is 
constant along particle trajectories: 

f(v,z,t) = f (v± (E,4)), 



(5) 



where zi. 



{z,t, E]U(z)) is the initial coordinate of 



a particle of total energy E which at time t is at the 
point z and has a velocity v±(E, z). Thus, z^ is a solu- 

-1/2 

dz . By defini- 



tionoft=(±) /* ± (E-U(z)) 
tion Sn R = f AvreB dv (/ - /,$ t=0) ) : 

f (v± (E,z±))-f (v ± (E,z)) 



5n 




V E - U{z) 



dE 



(6) 
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FIG. 1. Trapping potential. 

Here the sum is over particles moving to the right (+) 
and to the left (-), as in Fig. |l| The integration is over the 
resonant (negative) energies of trapped particles, U m < 
E < with U m being the amplitude of the potential. 

Let's first consider the short-time limit, r — > 0. Then 
the following approximations are valid. First (i), the 
effective resonance width § is AE res = (Av res ) 2 ~ 



l/(fer) 2 



oo, as r 



0. Second 



the particle ve- 



locity change is negligible \v(z, E) — v(zq, E)\ <C \/U(z), 
so that (iii) , the particle position is roughly proportional 
to time z^ ~ z ± t 2E j rrii. Finally (iv), the PDF re- 
sponse can be linearized (in the wave frame moving with 
v A ) as f (v) ~ f (v res ) +v fd(v rea ). Then Eq. @ may 
be estimated as 
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u(z±Ty/2E/r, 
U{z')dz'\ . (7) 



Here, we first used (i) to extend the integration over 
AE res to ±oo, and then (iii) and (i) to expand the de- 
nominator in (JsJ) for E ^> U. Finally, we took the r — > 
limit. Compared to Eq. (Q), the particle operator /C is 
replaced by the Hilbert operator, H, given by the ex- 
pression in curly brackets. It is nonlocal and satisfies the 
time-reversibility condition TtH = — 1. The effective dis- 
sipation coefficient is simply x\\ = 7r /o( 1 ' J 4)/ m i- Thus the 
KNLS equation (|l|) is recovered 

To treat the r — > oo limit, we recall that for the times 
t > tnl steady-state waveforms (discontinuities) have 
formed. Thus, particles are trapped in these adiabat- 
ically changing potentials. Hence, we may employ the 
virial theorem, which states that for any finite motion 
in a potential U(z) — U(z) — U m [i.e., U(z) > 0] the (pe- 
riod) averaged kinetic and potential energies are related 
by 



2(K(z))=n(U(z)). 



(8) 



Here U(z) is a homogeneous function of its argument of 
order n, i.e., U(az) = a n U(z). The resonance width is 
easily estimated to be Av tr — \j2\XJ m \jmi with \U m \ ~ 
B^/Stthq. Thus, for weak nonlinearity, the resonance is 
narrow. 



Av tr 

VA 



(9) 



Hence, an expansion of the PDF is valid, so that 
(in the wave frame): /o(u±) — fo(vA) ± w±/oGu) + 
(u±/2) fd'(vA). With this in hand and using Eq. (|8|) and 
(U) + (K) = E, we calculate the resonant particle con- 
tribution, Eq. (||): 



T — *00 

n 

n + 2 



(\U m \-\U(z)\) 



3(n + 2 



■\U(z)\ 



(10) 



Note that the term oc /q(wa) [ v ( z o ) ~ v ( z o )] vanishes 
identically because (U(zq )) = (U(zq )}. Thus damping 
is absent. Since (tC)(fC) KK — —1, we can, however, 
only estimate [from Eq. (0)] the coupling constant to be 
XII ~ /o("A)\/2/mj. The index n is formally not defined 
for an arbitrary potential. One may, however, estimate it 
comparing the calculated bounce period in the homoge- 
neous potential and "actual" one determined numerically 
for a known U , i.e., 



AE) = \E\ 



V E - U(z') 



(11a) 
(lib) 



It is interesting that the limit n —> oo encompasses two 
frequently encountered shapes of a wave packet, namely 
the solitonic and rectangular (i.e., deep narrow well) 
forms. In fact, for these cases as well as for any rather 
anharmonic potentials (n ^> 2) the resonant particle re- 
sponse ([lO]) is independent of n and takes on a very simple 
form: 



fa(vA)J-x\/W(*)\(\U m \-\U{z)\). 



(12) 



Thus, in the long-time limit, r 3> r tr , the damping rate 
vanishes due to phase mixing. Nevertheless, the resonant 
particles still contribute the wave dynamics, in that 



(5n R )~f '(v A )\bf 



(13) 



thus determining a new nonlinear wave equation. 

To estimate the number of trapped particles, we use 
a BGK-type (Bernstein-Green-Kruskal) approach [ pT| . 
This allows us to find the PDF such that the wave of a 
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FIG. 2. Asymptotic, r — > oo, particle distribution function. 

given profile is not dissipated by the Landau mechanism. 
In our problem there is only one resonance at v ~ va, 
since all modes are coherent. At large times, particle 
bounces result in flattening of the PDF at resonance so 
that Jq(va) — > 0. The height of this plateau (i.e., the 
phase density of trapped particles) depends on the wave 
evolution at earlier times. We take the "unperturbed" 
(trial) PDF as superimposed plateau and Maxwellian, as 
in Fig. |: 

fo(v) = F m (v) + [f p - F m (v)} Q VA (Av) , (14) 
where F m (v) is Maxwellian, the O-function is defined as 



Q VA (Av) 



1, if (v A - Av) <v<(v A + Av); 
0, otherwize; 



and f p is the constant to be determined. The coefficient 
f p has a simple meaning of the phase density of trapped 
particles after the plateau has been formed. Thus, the 
state with f p > F m (vA) corresponds to a clump on the 
PDF and that with f p < F m (vA) corresponds to a hole. 
The kinetic equation for a perturbation of the PDF is: 



(-iuj k + 7fc + ik\\v\\) = ik\\U kl 



dfo(v) 
dvn 



(15) 



By definition, Sn^^ = f /^fcdu. Then, for -y k <C LUk — 
k»VA and Av/va -C 1, we obtain: 



Stir = e 



k 



X 



-(ijk/k\\)F^v A ) + [F m (v A )-fp] 
~tl + kf.Av 2 



(16) 



Looking for the stationary solution, j). = 0, of the general 
KNLS equation (|j) and neglecting dispersion, we have 
d z [bSnNji + bSriR] = 0. Consequently, 



mib|b| s 



f P = F rn (v A ) + —2 

m 2 2vi 



(17) 



with Av = Av tr — va(B_i/B ). F m (vA) is the parti- 
cle phase density in the absence of trapping. Recalling 
that mi and m,2 are functions of [3 and x\\ °t f'o( v A) 
^ , we conclude that there must be an under-population 
of trapped particles [fo < F m (v A )] in a low-/? plasma 
((3 < 1) and an over-population [fo > F m (vA)] in a high- 
(3 plasma (/3 > 1). 

Finally, consider the there is weak wave damping not 
associated with Cerenkov resonance (e.g., as in ion- 
cyclotron or collisional damping). Then the wave am- 
plitude will slowly decrease, keeping resonant particles 
trapped. The following adiabatic invariant is thus con- 
served: 



J 



p\\dz ~ const, 



(18) 



i.e., (|w|||)(z2 — z\) ~ const. From Eq. ([!]), one can 
estimate Az ~ (fl c /vA){B±/ Bq)" 2 . Hence, Awy ~ 
(B±/Bq) 2 . The resonance width is, however, Avt r ~ 
(B±/B ). Thus, 



Av 



Av tr 



Bj_ 
Bo' 



(19) 



We thus obtain the trapped particle phase density: 



that is, the trapped particles will condense near the bot- 
tom of the potential well, as the wave amplitude de- 
creases. This results in a decrease in the effective index 
n, which approaches the asymptotic limit n — * 2. The 
BGK analysis given above is, however, then no longer 
applicable. It should be emphasized that trapped parti- 
cles condense in the bottom of the potential, rather than 
de-trap from it, as naively suggested in Ref. Jl2|] . Thus, 
no asymptotic, power-law damping exists in this case. 
Obviously, our considerations above are rather generic 
and valid for a wide class of nonlinear wave systems with 
quadratic nonlincarity and higher, and thus call the va- 
lidity of the results of Ref. |l^] into general question. 

To conclude, we have shown that the effects of the 
nonlinear PDF modification by a high-amplitude Alfven 
wave significantly modify the dynamics of such a wave. 
Even when phase mixing is efficient enough to quench 
linear Landau dissipation, trapped particles produce fi- 
nite a response which modifies the wave nonlinearity. 
The equation which explicitly describes the evolution 
of quasi- stationary Alfvenic discontinuities and asymp- 
totic (t — > oo) dynamics of nonlinear Alfven waves, Eqs. 
([| |l(j), has been obtained, this result constitutes the 
extension of the well established DNLS-KNLS theory of 
quasi-parallel nonlinear Alfven waves to the strongly non- 
linear regime of particle trapping. The phase density of 
trapped particles has been shown to be controlled by the 
value of plasma (3, as well as wave amplitude. 
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